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Abstract 



^ . We study a possible transition between symmetric nuclear matter and the diquark 

\ Bose-Einstein condensate (BEC) matter at zero temperature. We find that chiral 

' restoration transition is first order and coincides with deconfinement. We investigate 

\^ ■ various possible coexistence patterns which may emerge from the first order decon- 

\ finement phase transition by assuming different values for the critical deconfine- 



ment chemical potential. If deconfinement takes place at higher chemical potential, 
there exists a mixed phase of nuclear and chirally restored diquark BEC matter. 
This coexistence region extends over a large density region for a bigger diquark 
BEC or a weaker diquark-diquark interaction. For model parameters with heavy 
diquark in vaccum, phase transition to diquark matter becomes of second-order. 
P5 . We also show that in the case of precocious deconfinement, droplets of nucleons and 

droplets of chirally restored Bose-Einstein condensed diquarks coexist surrounded 
by non-trivial vacuum. We show that a early deconfinement and a weak repulsive 
diquark-diquark interaction soften the equation of state. We propose a scenario in 
5^ \ which nuclear matter saturates due to the formation of the diquark BEC and de- 

confinement phenomena. In this picture, instead of repulsive vector-meson exchange 
the compressibility of the equation of state is related to a repulsive diquark-diquark 
interaction. In general, we emphasize the importance of a diquark BEC phase at 
rather low density before quark BCS-pairing transition. 
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1 Introduction 



Exploring the phases of nuclear and quark matter has recently been subject 
of intensive investigation both theoretically and experimentally. The phase 
structure of QCD is relevant for a variety of phenomena, from heavy-ion colli- 
sions to neutron stars and the early universe. In particular, such studies may 
reveal the existence of matter in a new state which may be described by the 
elementary fields of QCD, i.e. quarks and gluons. 

One of the most interesting question in QCD phase diagram is to find the 
phase boundary between the hadronic and quark phases. It has been well do- 
cumented that at zero baryonic density, chiral restoration and deconfinement 
phase transition occur at the same temperature [1] . Our current understanding 
of QCD phase structure at finite density is still rudimentary. Ab initio lattice 
Monte-Carlo simulation cannot be directly used due to the complexity of the 
sampling weight, the so-called sign problem. On the other hand, there is no 
alternative approach to accommodate all relevant non-perturbative features 
of QCD and describe both phases in a unified way. Various lattice simulation 
techniques have recently been proposed in order to circumvent the lattice sign 
problem (for a recent review see Ref . [2] ) . But all of them still suffer from sys- 
tematic uncertainties which made them reliable at very low chemical potential. 
At low temperature and large densities, a number of model calculations appear 
to agree on a first order phase transition between nuclear and quark matter 
[3,4]. Therefore, one expects that a first order phase transition at finite density 
will end in a second order endpoint and then turns into a rapid crossover at 
small chemical potential [2]. If confirmed, this indicates that as we approach 
to lower temperature and finite baryonic density, the phase transition becomes 
stronger. 

At the nuclear matter saturation point the Fermi momentum kp and the pion 
mass TTiTr are of comparable scale kp ^ 2171^^ ^ 260 MeV. This implies that at 
the densities of interest in nuclear physics po = 0.15 fm^'^ = 0.45 m^, pions and 
chiral dynamics must be important. The methods of chiral perturbation theory 
exploit this observation and have been successful to describe the nuclear mat- 
ter properties [5] . By estimating the volume occupied by a single nucleon from 
its mean square charge radius rrms = \f(j^, one may conclude that baryonic 
matter is dilute up to about 3po, since it can satisfy the diluteness condition 
'"rmsP^^"^ <^ 1. At higher density nucleons start to overlap. Nuclear matter is 
dilute and at the same time it is strongly correlated since ap^/^ ^ 1, where a is 
the scattering length. On the other hand, since the Fermi momentum exceeds 
the QCD scale parameter Aqcd ~ 200 MeV, multiple scattering between nu- 
cleons on the Fermi surface start to probe distance of the order ^ 1 fm and 
the nucleonic substructure becomes visible. The Fermi momentum at satura- 
tion point is in fact of the order of kp ^ Aqcd- The confinement scale can 
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even become one order of magnitude smaller at moderate baryonic density 
[6]. Therefore, it is legitimate to ask whether confinement and nuclear matter 
saturation mechanisms are interconnected. The delta-nucleon mass splitting 
A Ki 290 MeV is also comparable to the Fermi momentum kp ai nuclear mat- 
ter saturation. It is known that the main mechanism behind the delta-nucleon 
mass splitting is related to their different internal diquark configurations [8]. 
Namely, the delta is made of only axial vector diquarks while nucleons can 
be constructed from both the scalar and axial vector diquarks. This indicates 
that perhaps diquark dynamics at nuclear matter saturation point might be 
also relevant. The relevance of diquarks as an efficient way to incorporate the 
non-perturbative feature of QCD and to simplify the underlying dynamics 
has been well-known for years (cf. [7,8] and references therein). For example, 
baryons properties in vacuum have been successfully described in the covariant 
diquark-quark picture [8] . Recently, it has also been shown that nuclear matter 
properties can be described by simulating the nucleons as relativistic bound 
states of the diquarks and quarks [9,10]. At high density diquark concept is 
also crucial to describe the color superconductivity phenomenon [11,12]. On 
the other hand, if diquarks as bound states exist, they can in principle un- 
dergo Bose-Einstein condensation since they are bosons. Therefore, it seems 
to be natural to assume that baryonic matter at some density undergoes phase 
transition to diquark BEG matter. One of the motivations behind this paper 
is to explore such a transition. 

Based on the above scaling argument, it seems that the most important non- 
perturbative features of QCD might be involved in the nuclear matter satura- 
tion density. In this paper we attempt to address such a possible connection 
between the nuclear matter saturation properties and chiral restoration and 
the deconfinement effect at zero temperature. We simulate the quark decon- 
finement at finite density by allowing that the active degrees of freedom change 
from nucleons to diquarks at high density. This may lead to a first order 
deconfinement phase transition, in accordance with indications from earlier 
studies [3,4]. We calculate the nuclear and diquark equations of state (EoS) at 
zero temperature. In particular, we will investigate a possible diquark Bose- 
Einstein condensate at about the nuclear matter density. Our result leads to a 
non-convential scenario for the stability of the nuclear matter. We show that 
the nuclear matter is energetically favored to undergo the diquark BEG phase. 
We study the emergence of various coexistence patterns by allowing the de- 
confinement to occur at different critical baryonic chemical potential. We will 
show that it is possible that: all three transitions chiral restoration, the de- 
confinement and the liquid-gas transition are coincident at zero-temperature. 

This paper is organized as follows: In section 2, we first study the diquark 
BEG dynamics. In section 3, we introduce the hadronic part of our model 
and calculate the nuclear and diquark EoS. In sections 4 and 5 we investigate 
the transition between diquark BEG phase and nuclear matter phase at high 



3 



and low baryonic density, respectively. Some concluding remarks are given in 
section 6. 



2 Diquark BEC and spectrum of Goldstone modes 



With any weak attractive interaction fermions exhibit superfluidity charac- 
terized by an energy gap from the BCS pairing mechanism. (For reviews on 
color superconductivity, see Refs. [13,14]). However, if the attractive interac- 
tion is strong, the fermions can form bound bosonic molecules which condense 
into the bosonic zero-mode at some critical temperature. This is called Bose- 
Einstcin condensation (BEC). These two different mechanisms are in fact two 
different physical manifestations with the same order parameter and are ex- 
pected to be connected by a crossover ^ . In the weak-coupling region at high 
density the correlation length between two quarks is bigger than the interquark 
distance or the pair size. However, one may expect that in strong-coupling re- 
gion at rather low density due to large pair fluctuations, the correlation length 
becomes smaller than the pair-size and quark matter undergoes the diquark 
BEC phase [16]. Therefore, at the density around the deconflnement transi- 
tion, the diquark bound states may be the relevant degrees of freedom. Here, 
we are interested in the physics of low density (strong-coupling region) where 
the BEC phenomenon is more relevant. Therefore, wc assTimc that the di- 
quarks bound states are already formed and we concentrate on Bosc- Einstein 
condensation. The diquark matter has been already studied in Ref. [17] but 
without taking into account the diquark BEC effect. We use the following 
Lagrangian density for a diquark fleld 

= (^0 + i//)$^(ao - i//)$ - (a,$) - m^$t$ _ g (^^^^y + . . . , (i) 



where ii is the chemical potential with respect to the U{1) charge and 
denotes the diquark mass. The coupling g gives the repulsive diquark-diquark 
interaction. We assume that the higher-order scattering terms can be ignored. 
$ is a color 3 complex scalar fleld which represents the diquark degree of 
freedom. 







1 

71 



^21 



'12 



'22 



031 + «032 



(2) 



where the complex diquark flelds are parameterized in terms of a possible 
condensate 0o — (0) (due to symmetry, the direction of condensate is at our 



^ This phenomenon has been experimentally examined in various ultracold gases 
of fermionic alkali atoms using the Feshbach resonance techniques [15]. 
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disposal) and fluctuating flclds (pai ^-nd 0q2 with a = 1,2,3. We have chosen 
the diquark condensate to be in a = 3 direction. The Lagrangian Eq. (1) 
possesses a SU{3) x U{1) global symmetry which spontaneously breaks down 
to SU{2) X U{1)q in the presence of the condensate. This is exactly the same 
symmetry pattern which underlies the 2SC phase [13,14]. Having substituted 
the diquark field parameterization Eq. (2) into Lagrangian Eq. (1), one can 
represent the diquark Lagrangian Cd Eq. (1) as a sum of the condensate La- 
grangian Cc, free piece Cp with effective masses mi and m2 and the effective 
interaction part £/ among various modes (j)ai and (f)a2, 



c '■ 



= ^031 {d^^d, - ml) 031 + ^011 {d^d, - ml) cj>n 
Ihi {d'-d, - ml) 021 + E (^'^«2 (d'-d, - 



+ 2 



ml) 0„2 



^«i 



dt 



dt 



- //^ + 2g(f)l + gY^ (0al0al + 0a20a2) 



- I ('/'al^al + 0a20a2)^ , 



000: 



'31 



(3) 



with 



ml — —(J? + + 6(70^ 



0' 

mi = -//^ + + 2^02. (4) 

The inverse of free diquark propagator corresponding to Lagrangian Eq. (3) 
is a 2 X 2 matrix in (0ai, 0a2) space for different color directions, and in 
momentum space it reads 



iV^^{uj,p)^iV^^{u,p) = 



2 2 2 

UJ — p — m2 



iV^^{(jj,p) 



2 2 2 

o;"^ — — mf 



■2iiMjj uj"^ — p^ — ml j 
2iiiu! 



\ 



-2iiiu! cu^ — p^ — m^ 



(5) 



2/ 



The excitation spectrum is found from the diagonization of the two-point 
Green function or, equivalently, by solving the conditions Det(P~li_2,3) = 0- 
One finds that six degrees of freedom of the complex anti-triplet diquark field 
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split into a degenerate pair of eigenvalues for a = 1,2 and a singlet pair 0i>3± 
for {a = 3) (we only consider positive eigenvalues, the negative ones can be 
treated similarly). 



a;i,2± = iyw + V + + "^1 (a = 1,2), 

ujs± = ^p^ + ^ (mf + m|) + 2//2 ± ^\/(mf + m| + 4//^)^ + 4:{Aii'^p^ - m\ml) 

(a = 3). (6) 

It is easy to show that in the absence of a condensate 0o = the splitting 
of the three color directions goes away, as it should. For = and no sym- 
metry breaking > excitation energies are all identical uj — ^Jk^'+m^ as 
expected. 

The classical potential has a trivial solution 0o = and a nontrivial minimum 



= (7) 

In the case that < 0, symmetry breaking would even occur in vacuum since 
//^ — is then greater than zero. We assume, however, that > which 
means no symmetry breaking at zero baryon density. Only when //^ > we 
have condensation and the trivial 0o = is no longer a minimum of potential 
energy. The system develops an instability with respect to formation of a 
diquark condensate. In this ground state, the initial SU{3) x U{1) symmetry 
spontaneously breaks down to SU{2) x U{1). Using Eq. (7), the effective masses 
in Eq. (4) simplify to mf = 2(/i^— m^) and m2 = 0, and the dispersion relations 
in Eq. (6) become 



UJ3± = + 3^2 _ ^2 ^ ^ (3^2 _ ^2)2 ^ ^^2^2 _ (g) 



The small momentum expansion of the dispersion relations in far infrared 
region is given by 
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"^1,2+ 



<^1,2- 




(10) 



(9) 



(11) 




(12) 



Prom the above relations it is noted that the excitation 001^2+ and have 



gaps identical to 2yU and yGyU^ — 2m^, respectively. The Nambu-Goldstone 
(NG) bosons can be identified with excitations which are gapless as their mo- 
mentum goes to zero. Therefore, there are only 3 gapless NG bosons in the 
spectrum. Due to symmetry consideration one may expect 5 NG bosons. If the 
conventional counting of NG bosons would apply, the number of NG bosons N 
is equal to the number of broken generators. However, as has been previously 
noticed in Lorentz-noninvariant systems the number of NG-modcs can be less 
than the number of broken generators [18]. A more accurate counting role is 
due to the Nielsen-Chadha theorem [19] which states that if rii and 712 are 
the number of gapless excitations with infrared dispersion laws a;i p and 
UJ2 — p^, respectively, then the number of NG bosons N < ni + 2n2. Therefore, 
there is no contradiction since two of three gapless modes Eqs. (10,12) have 
quadratic dispersion laws for long wavelengths. It is interesting to note that 

2 

at the BEG phase transition point when /i rud we have uj^- — f^, and 
all three gapless excitations become degenerate with the quadratic dispersion 
law. Therefore, based on the naive continuity of spectrum on the boundary of 
phase transition, one might have expected less NG bosons. One should also 
note that since we have gapless mode with quadratic dispersion laws rather 
linear, the Landau criterion [20] for superfluidity fails ^. Dissipationless flow 
in homogeneous bosonic system is one of the manifestation of the interplay- 
ing microscopic interaction and the Bose-Einstcin condensation. The question 
whether a Bose-Einstein condensed diquark gas in a more realistic model re- 
ally displays the Landau criterion for superfluidity remains to be answered 
elsewhere. 

The thermo dynamical potential ild(/i, m^) at diquark chemical potential /i is 



According to the Landau criterion, a fluid moving with group velocity v is su- 
perfluid if v < Vc = Min C^ ) where the minimum is taken with respect to all 
excitation modes and the momentum p. If system has a quadratic dispersion of the 
low-energy excitation modes, then Vc = 0, and Landau criterion can not be satisfied. 
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obtained through the partition function Z 

Q,(/i,m,) = P^t^gi/ci^./:.^ (13) 

where V denotes the spacetime volume of the system. The one-loop contribu- 
tion to the effective potential evaluated at the classical minimum Eq. (7) is 
obtained as 



r dio 




1 2^ J 


(27r)3 y 



2 In Det iX)i-^(a;, p) + In Det iP3"^(cu, p) 



4y 



(2 



(14) 



where Vil is defined in Eq. (5) and evaluated at the minimum of the classical 
action. The a;i,3± is defined in Eq. (8). In the second line in the above equation, 
wc have performed the integration over energy using contour integration. The 
first term 0,^^^ in Eq. (14) is the mean-field thermodynamical potential evalu- 
ated at the minimum of the classical action. The momentum loop-integrals in 
the diquark effective potential Eq. (14) are divergent and have to be regular- 
ized. For simplicity we use a sharp cutoff A/j in three-dimensional momentum 
space. Note that scalar field theory given in Lagrangian Eq. (1) is renormal- 
izable (when higher order interaction terms are ignored), therefore any cutoff 
dependence can be removed by some kind of renormalization scheme. Here, 
we assume that the model has a natural ultraviolet cutoff which is defined for 
the vahdity of the model at high-energy. We will investigate the imphcation 
of various choices for the cutoff. 



One can show that the zeroth and first-order derivatives of the thermody- 
namical potential Eq. (14) with respect to chemical potential fj, at fj, — rrid 
are continuous, but the second-order derivative at /i — rud is discontinuous. 
This indicates a second-order phase transition between normal and diquark 
condensate phases. The charge density p and the speed of sound Vg can be 
obtained from the thermodynamical potential 



(16) 

Using the above relation, the speed of sound at mean-field approximation is 
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obtained from the first term in the thermodynamical potential Eq. (14), 




(17) 



Note also that this is identical to the coefficient of u^- in Eq. (12). At very 
high density fi ^ we recover the standard relativistic result "^s = 
independent of the model parameters. This indicates that at very high density 
the NG boson propagates as in a relativistic medium. 



3 DiquEirk BEC and the nuclear matter equation of state 

We incorporate the deconfinement effect by allowing that the active degrees 
of freedom change from baryons to diquarks at a critical baryonic critical 
chemical potential /i^. For description of the hadronic phase, we use the chiral 
a — CO model with exact global SU (2) x SU (2) symmetry, which contains a 
pseudoscalar coupling between pions and massive isoscalar vector field uj^ to 
nucleons. 



where the field tensor is defined = d^uOij—dyUJ^ and ^5 denotes the baryonic 
chemical potential. The nuclcon mass Mjv and the vector-meson mass at 
rest are generated through spontaneous symmetry breaking by the cr-field 



The energy density for the hadronic phase can be obtained from Eq. (19) in the 
conventional mean-field approximation. Having used the equation of motion 
for the field, we obtain the energy density as 




(19) 



m„((7) = g^a. 



(20) 



Sh = V{<t) + 



9vnPb 




2ml{a) 



(21) 




(22) 
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where we defined ^^{k, a) — ^k^ + M^{a) and a ( and ctq) denotes the mean- 
value of the sigma field in the nuclear matter medium (and in the vacuum). We 
have also included the corresponding Dirac-sea of the nucleon up to a ultravi- 
olet cutoff Atv- In principle, in renormalizable models such as the linear sigma 
model, divergences can be absorbed into the coefficient of interaction terms 
by employing a renormalization scheme. However, it has been shown that for 
this kind of models with heavy fields, such terms give rise to unnaturalness ^ 
indicating that quantum vaccum is not adequately described by long-range 
degrees of freedom [21]. We have recently shown that nuclear matter proper- 
ties can be better reproduced if one treats the vacuum loop explicitly [24] . We 
discard the short-range physics which cannot be described by long-range de- 
grees of freedom by the nucleonic ultraviolet cutoff Ajv- The baryonic density 
Pb is related to the Fermi momentum kp hy 

At any density pb, £h is stationary with respect to the mean- value of the scalar 
field. Therefore, the scalar mean-field a is obtained via the self-consistency 
conditions: 

^ = 0. (24) 
dcj 

This equation is numerically solved for every point of density. The baryonic 
chemical potential ph can be obtained through the Hugenholtz-van Hove the- 
orem via the energy density 

/^;'-$^ = 4rT + M^i^,^), (25) 
dp m^{a) 

where the value of a for a given density is obtained via Eq. (24). Then the 
thcrmodynamical potential Oh and the hadronic pressure are obtained 
through a Legendre transformation and can be written as 



2ml{a) 



dk k {ujN{k,a) — ujN{k,ao)) ■ (26) 



^ Furnstahl et. al. [21] have shown that the one-baryon-loop vacuum contribution 
in renormahzcd models like the linear sigma model and the Walecka model gives rise 
to large unnatural coefficients based on the "naive dimensional analysis" proposed 
by Georgi and Manohar [22]. 
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For the deconfincd phase we consider matter made of diquarks coupled to a 
and TT and we use the diquark Lagrangian introduced in the previous section. 
The new input here is that we couple the diquark field to the scalar field in 
order to generate the diquark mass dynamically via spontaneous symmetry 
breaking. 

£q - {^^+^^x)^\^o-^^x)^-{^i^^) {di^)-{gd'jf^^^-g ($t$)'+£^^. (27) 

The Lagrangian vC^.^ is defined in Eq. (19). We incorporate the diquark BEC 

efi^ect in the fashion presented in the previous section, namely by employing 
the diquark field parametrization Eq. (2). We treat the diquark fields at the 
one-loop level. Therefore, the thermodynamical potential ilq and deconfined 
pressure Pq are given by 

Qq = -Pq = Qd(/x, ma{a)) + V{a) - ^ j^" k^dk^k^ + ml{ao), (28) 

where the diquark mass is related to the cr-field: 

rridicr) = QdCr. (29) 

The fid and V{a) in Eq. (28) are defined in Eqs.(14) and (22), respectively. 
The last term in Eq. (28) is the constant diquark contribution at zero density 
which shifts the minimum of the vacuum thermodynamical potential to zero. 
Note that the diquark mass now dynamically changes in medium through the 
in-medium cr-field. This does not alter the possible formation of the diquark 
BEC as far as /^^ > where ii is the diquark chemical potential. In contrast 
to the Bag model and similarly to the NJL model, our diquark model dy- 
namically generates a density-dependent bag constant and density dependent 
effective constituent diquark masses. With the diquark thermodynamical po- 
tential Eq. (28) and the density Eq. (15) one can construct the energy density 
through the Legendre transformation. 



1 r^D p^dp 

+ j_ p^dp ( (p' - ^li^)) ^3+^3- + - PV ^ 

- ^ j^'' edk^k' + mjicTo), (30) 

where for the right hand side we have made use of the definition of p(/i) and 
Qq given in Eqs. (15,28). The uj-^-j. and V^a) are defined in Eqs. (8) and (22), 
respectively. Using Eq. (8) we obtain the following relations which facilitate 
the computation 
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+ = \/2^p2 + 3/^2 - ml{a) + p^J + 2{iJ? - ml{a)) . (31) 

The in-medium mean- value of the scalar field and the corresponding diquark 
mass is obtained similarly to the hadronic sector by minimization 

^ = 0. (32) 
d(7 



The diquark chemical potential and density (//,p) are related to baryonic ones 
(a*6, Ph) by: PL — l/ift and p — jpi,. The binding energy per baryon is defined 
as 

^ = ^ - M^(ao) = + p,- M^(ao), (33) 

A Pb Pb 

where eh,q and ph,q are defined in Eqs. (21,30) and Eqs. (26,28), respectively. 



4 Nuclear-diquark phase transition at high density 



We incorporate the deconfinement effect by demanding that as we increase the 
density in the hadronic phase at a critical Fermi momentum kp corresponding 
to the baryonic chemical potential f^l, the baryonic degrees of freedom are 
replaced by diquark ones while keeping the description of the auxiliary scalar 
fields intact above the transition. We do not couple the vector-meson field to 
diquarks. The exact value of the deconfinement chemical potential fil is not 
known. Therefore, we will consider the implication of different choices for pi 
and kp. 

In the hadronic phase, our model contains six unknown parameters A^r, m^, 
A, Qn, Qvn and g^. We take the hadronic ultraviolet cutoff Ajv = 324 MeV [24]. 
In this model the acceptable range of the ultraviolet cutoff Ajv is rather low 
[24]. This empirical fact has been also observed by other authors within the 
NJL model with nucleonic degrees of freedom [25] . This is in contrast with the 
case that the model is defined in terms of quark degrees of freedom where the 
allowed range of the cutoff is typically larger. Due to spontaneous symmetry 
breaking the scalar mean field acquires a nonzero vacuum expectation value 
which is equal to the pion decay. The pions are then the massless Nambu- 
Goldstone bosons. For the given A;v = 324 MeV, one can determine the scalar 
potential couplings and A so that the mean- value of the scalar field becomes 
identical to the empirical pion decay at vacuum ctq = 93 MeV. We take the 
empirical nucleon mass in vacuum Mjv(o'o) = 939 MeV and vector-meson 
mass in vacuum m^(cTo) = 783 MeV. These choices fix the coupling Qn ~ 
-^at/Ztt = 10.1 and — m^j — 8.42. The vector couphng g^^ is treated 
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parameter Aiv(MeV) mg(GeV^) A gyn g-n Ov 

set NM 324 +0.026 25 6.54 10.1 8.42 

Table 1 ' 

The value of various couplings in hadronic phase for set NM is given. This parameter 

set reproduces the empirical saturation point {Eb/A = —15.75 MeV, po = 0-15 fm^). 

as a free parameter and is adjusted in the medium in such a way that the 
EoS reproduces nuclear-matter saturation properties, namely a binding energy 
per nucleon Eb/A = —15.75 MeV at a density corresponding to a Fermi 
momentum oi kp = 256 MeV = 1.3 fm~^ [23]. In this way all parameters in 
the confined phase are fixed uniquely (see set NM in table 1). 

The sigma mass can be obtained from 

- (34) 

where the derivatives are evaluated at the scalar mean-field a solution of the 
gap equation. The corresponding sigma meson mass nia in vaccum for the 
above parameters (set NM in table 1) is 810 MeV. In order to estimate the 
stiffness of the nuclear matter EoS, one may compute the compression modulus 
at the saturation density. It is defined as 

K = -9-^ = Qpb^T- (35) 
dpb dpi 



The compression modulus at saturation density for the parameter set NM is 
K — 455 MeV. Note that in the original a—u model of Boguta the compression 
modulus of nuclear matter is unacceptably large about K = 650 MeV [26]. 
Here, due to the inclusion of the Dirac-sea we are relatively able to reduce 
the compression modulus (see also Ref. [24]). This value is still bigger than 
the empirical one which is beheved to be about K — 200 — 300 MeV [28,29]. 
Nevertheless, here we use only the nuclear matter compression modulus as a 
relative measure of stiffness of the deconfined EoS. It has been very difficult to 
overcome the unphysically high values of compressibility in the a—uj model and 
in the Walecka model. In the next section we will introduce a non-conventional 
scenario in which the nuclear matter compressibility can be lowered. 

For zero temperature, two phases of hadrons and diquarks can be in equilib- 
rium if their chemical potentials and pressures are equal (Gibbs' conditions): 

PM)=Pn{pX), (36) 

where Ph and Pg are defined in Eqs. (26) and (28), respectively. First, we as- 
sume that the transition to the diquark BEG matter occurs at the baryonic 
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Pnrfnnpfpr ( ii* — 1914 5 MpV^ 


qpt Al 

t JT. J. 


<5pt A 9 




opt A A 


a 


14.1 


13.47 


12.73 


11.23 


9d 


6.85 


6.7 


6.5 


6.0 


md(p = 0)(MeV) 


637.1 


623.1 


604.5 


558.0 




3.0 


3.0 


3.0 


3.0 


Pb/PO 


3.0 


4.0 


5.4 


8.6 



Table 2 

The parameters and 5 in the deconfined phase for the sets A1-A4 for the critical 
deconfinement chemical potential /x^ = 1214.5 MeV. We assume as ultraviolet cutoff 
for all parameter sets Ad = 900 MeV. The baryonic density at the transition point 
in hadronic and diquark parts in units of the normal nuclear matter density 
po = 0.15 fm^ are also given. 

chemical potential — 1214.5 MeV which corresponds to the baryonic den- 
sity = 3po from the hadronic side. For the deconfined phase, we have 
three unknown parameters: the ultraviolet cutoff A^, the diquark-diquark 
coupling g and the diquark-scalar field coupling gd- We assume the decon- 
fined ultraviolet cutoff to be Ad = 900 MeV and fix the other two unknown 
parameters by Gibbs' conditions Eq. (36). We require that at the phase tran- 
sition the baryonic chemical potential of both hadronic and diquark parts to 
be equal ^5 = |/x = /i^ = 1214.5 MeV. Having obtained the nuclear mat- 
ter EoS, one can immediately read off the hadronic pressure at the critical 
deconfinement chemical potential: Ph{l^X — 1214.5 MeV) = 89.8 MeV/fm^. 
Therefore, the Gibb's condition for equilibrium of two phases is satisfied if 
PhifJ-b) = Pqif^b) = 89-8 MeV/fm^. This gives a non-trivial constraint on the 
diquark-model parameters. There exist many solutions which satisfy the above 
requirements, each of them corresponds to a different BEC condensate. The 
values of parameters g and g^ for four sets A1-A4 with a fixed Ajj = 900 MeV 
are given in table 2. We have also given in table 2 the value of the diquark 
mass in vacuum using Eq. (29). Note that the values of the scalar diquark 
mass in vacuum for various parameter sets given in table 2 is within the ac- 
ceptable range which has been estimated by many authors in order to describe 
the baryon properties in vacuum [8,7,27]. 

In Fig. 1, we show the mean- value of the scalar field a with respect to the 
baryonic chemical potential fib- The solid line denotes the mean-value of the 
scalar field in the nuclear matter medium when the deconfinement to diquark 
matter is not incorporated. For all diquark parameter sets (given in table 
2) chiral phase transition takes place at the critical deconfinement chemical 
potential = 1214.5 MeV. We will show later that this will also be the case 
if the critical deconfinement chemical potential is taken lower value. 

In Fig. 2 (right pannel), we show the diquark pressure Pq for sets A1-A4 (in 
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Fig. 1. In-medium pion decay constant as a function of baryonic chemical potential 
/Xfc. We use the parameter set NM given in table 1 for nuclear matter and parameter 
sets Al — AA given in table 2 for diquark matter. 

table 2) and the hadronic pressure Ph for set NM (in table 1) as a function 
of the baryonic chemical potential /ib. For comparison we also display the 
nuclear matter solution for the parameter set NM when the transition to a 
diquark BEC is not incorporated. It is obvious that the EoS becomes consider- 
ably softer if deconfinement is allowed. Therefore, the system is energetically 
favored to undergo phase transition. By comparing different diquark EoS ob- 
tained by using various parameter sets given in table 2, it is observed that a 
bigger diquark-diquark coupling g leads to a stiffer EoS, whilst a bigger cou- 
pling g leads to a smaller diquark condensate, see Eq. (7). Therefore, the bigger 
the diquark BEC condensate the softer the EoS will be. The inserted plot in 
Fig. 2 shows the first-order liquid-gas phase transition of nuclear matter which 
takes place at ix^ = 923.25 MeV. The first-order phase transition is manifested 
by the appearance of several branches of Ph{nb)- At a fixed baryonic chemical 
potential, only the highest pressure corresponds to a stable phase. The slope 
of curves P{fib) exhibits a jump in density. In Fig. 2 (left pannel), we show the 
baryonic density pb in units of nuclear matter density po as a function of the 
baryonic chemical potential, we have also plotted the critical deconfinement 
chemical potential fi^. In general, there are two separated density jumps with 
respect to the chemical potential: one corresponds to the standard nuclear 
matter liquid-gas phase transition and the second one corresponds to the de- 
confinement phase transition. The value of baryonic density below and above 
deconfinement is given in table 2. Depending on model parameters the jump 
at the deconfinement transition can change. The softest EoS set A4 results in 
a high baryonic density above the phase transition since the diquark BEC con- 
densate is bigger and consequently it can accommodate more baryon density. 
Note also that the set A4 has the lowest diquark mass in vaccum, see table 
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Fig. 2. Baryon number densities (left) in units of nuclear matter density 
po = 0.15 fm~^ and pressure (right) as functions of the baryonic chemical potential 
Hb for various parameter sets given in tables 1 and 2. The critical deconfinement line 
fj,^ = 1214.5 MeV is also shown. The insert shows the region of the nuclear matter 
liquid-gas phase transition. 

2. For the model parameter (set Al) with rather heavy di quark in vaccum, 
transition to diquark matter takes place without a density jump, i. e. the first 
derivative of the thermodynamical potential with respect to the chemical po- 
tential is continuous. This indicates a second-order phase transition between 
hadronic and diquark matter. In this case, although deconfinement and chiral 
phase transition coincides, their phase orders are different. 



In order to consider the possible coexistence regions which emerge from the 
first-order phase transition, we perform Maxwell constructions. In Fig. 3, we 
show the energy per baryon eh,q/ pb as a function of inverse baryonic density 
(Pb/po)^^ for the parameter set A4 given in table 2. It exhibits a typical 
Van der Waals-like diagram. Using the following thermodynamical consistency 
equation and Gibb's conditions, 

A (—) = -Pn,,. (37) 
dpb \ Pb J 

one can identify the coexistence region as a function of inverse density. There 
exist two separated coexistence regions which are denoted in Fig. 3 as coex- 
istence region 1 and 2. At low density < p^, < po iii the hadronic phase we 
have the typical coexistence between nucleonic gas and nuclear droplet. At 
rather high density 3po < Pb ^ pt (J^^ values of pf see table 2), there exist a 
coexistence region between nuclear matter and chirally restored diquark BEG 
matter. At zero temperature, the coexistence region between the hadronic and 
diquark BEG phase seems to be quite large in the density axis for a weaker 
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Fig. 3. Energy per baryon number as functions of inverse baryonic density in units of 
the nuclear matter density, i. e. {pb/po)~^ when the deconfinement chemical poten- 
tial is taken fi^ = 1214.5 MeV. The dotted lines represent the Maxwell-constructions 
which determine the regions of coexistence. The coexistence region 1 contains the 
nucleonic gas and nuclear droplets and the coexistence region 2 contains mixed 
nuclear and diquark BEC matter. 

diquark-diquark interaction or a bigger diquark BEC (see parameter set A4). 
Forcrand et al. [32] have recently found some indication for the coexistence 
region of the hadron and the plasma phases based on the canonical approach 
of lattice QCD (for 4-flavour and pion mass k, 350 MeV). Their extrap- 
olation indicates that the coexistence region might be quite wide in density 
0.50(5) B/fm^ ^ P ^ 1.8(3) B/ fm^. It is of great interest to see if this coex- 
istence region will survive in the chiral limit using a fully realistic simulation 
of QCD. 



5 Nuclear-diquark phase transition at low density 

In this section, we consider the possibility that the deconfining transition 
happens at a lower baryonic chemical potential. In particular, we study the 
extreme case where the deconfinement phase transition takes place right at the 
nuclear matter saturation Fermi momentum, namely at k*p = 256 MeV which 
is approximately equal to the QCD confinement scale Aqcd- This corresponds 
to a baryonic chemical potential /i^ = M7v(o"o) — Eb/A = 923.25 MeV. At 
the saturation point, the binding energy per nucleon is taken Eb/A = 15.75 
MeV. Above the deconfinement phase transition the vector-meson does not 
appear in the deconfined Lagrangian Cg Eq. (27) and one may be afraid that 
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Table 3 

The parameters ga and g in the deconfined phase for the sets B1-B4 for the critical 
deconfinement chemical potential fi^ = 923.25 MeV. The diquark ultraviolet cutoff 
for all parameter sets is taken Ad = 650 MeV. 

nuclear matter in this model would collapse [26]. However, we will show that 
the repulsive diquark interaction induces enough repulsion to reproduce stable 
matter similar to conventional nuclear matter. 

For the case of early deconfinement, a lower diquark ultraviolet cutoff Ao 
is allowed since we are interested in physics of rather low baryonic chemical 
potential. We assume as diquark ultraviolet cutoff A = 650 MeV. Later we 
will investigate the effect of choosing a higher value for the ultraviolet cutoff. 
The acceptable range of the diquark-scalar field coupling ga and the diquark- 
diquark coupling g is determined via the Gibb's condition. Since from the 
hadronic part at the saturation point wc have that the hadronic pressure 
Phil^h) = 0' require that also the diquark pressure Pq{nX — 923.25MeV) = 
0. In table 3 we give the values of the couphngs in four parameter sets B1-B4 
obtained by the above-mentioned procedure. 

In Fig. 4, we show the in-medium pion decay constant /^r = ex as a function 
of the baryonic chemical potential (left pannel) and baryonic density Pb/po 
(right pannel). The chiral restoration for all parameter sets given in table 
3 is first order and takes place right at the critical deconfinement chemical 
potential /i^. The solid line denotes the mean-value of the scalar field in the 
nuclear matter medium without deconfinement effect. It is well known that in 
the (7—uj model with dynamical 6i;-mass, chiral restoration does not occur since 
a tends to grow at high density. This can also be seen in Figs. 1 and 4, where a 
at higher chemical potential and density in nuclear medium (solid line) bends 
upward. In fact this feature of the model is closely connected to the stability 
of the system. The sigma field plays the role of the chiral partner of the pion 
and at the same time it is the mediator of medium-range nucleon-nucleon 
attraction and via the a;-mass also the short range repulsive is related to the 
sigma field. Therefore, as its mass becomes smaller the attraction between 
nucleons becomes stronger which may destroy the stability of nuclear matter. 
This effect is not present in non-chiral type models such as the Walecka model 
[23]. However, we will show that in our approach the chiral restoration does 
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Fig. 4. In-medium pion decay constant a as functions of density Pb/po (left) in unit 
of nuclear matter density po = 0.15/m~'^ and baryonic chemical potential pb (right). 
We have used parameter sets NM and set B2 given in tables 1 and 3. 
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Fig. 5. Thermodynamic potential fig (and fi/j) as a function of auxiliary variable 
(T (pion decay) for parameter set NM (for pb = 0) and for other values of pb the 
parameter set B2 is used. 



not obstruct the stability of the system. 



The thermodynamic potentials fig and as a function of the auxiliary vari- 
able a are shown in Fig. 5. In vacuum the minimum of the thermodynamic 
potential flh is located at the pion decay ex = 93 MeV. As we increase the 
chemical potential right at the critical deconfinement chemical potential fil 
where we switch to the diquark thermodynamic potential Qg, there is a strong 
first order phase transition from the vacuum to the chirally restored phase with 
a = 0. Once the phase transition has taken place, the minimum of f2g at a = 
remains with increasing chemical potential, but the lowest Qq configuration 
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Fig. 6. Baryon number densities (left) in units of nuclear matter density 
Pq = 0.15/m~'^ and pressure (right) as functions of the baryonic chemical poten- 
tial ph for various parameter sets given in tables 1 and 3. For comparison we have 
also plotted the nuclear matter solution (solid line) for the parameter set NM when 
transition to diquark matter is not taken into account. 

becomes more and more negative fig < 0, indicating an increase in pressure 
Pq > with density. At the critical chemical potential ix\ = 923.25 MeV, both 
minima of Qh and flq have zero pressure. One represents the non-perturbative 
broken vacuum at density zero and the other represents the chirally restored 
diquark BEC phase at non-zero density ph- Note that we construct the de- 
confinement phase transition in such a way that the chiral phase transition 
results. The chiral restoration is first order and coincides with the deconfine- 
ment critical chemical potential. 

In Fig. 6 (right pannel), we show the nuclear matter pressure Ph for set NM 
(in table 1) and the diquark matter pressure Pq for set B1-B4 (in table 3) as a 
function of the baryonic chemical potential. It is seen that above the critical 
deconfinement chemical potential the diquark EoS becomes generally softer 
than nuclear matter solution. As in the previous case (in Fig. 2), the behavior 
of the diquark EoS depends on the diquark-diquark coupling g and conse- 
quently on the strength of the diquark Bose-Einstein condensate. A bigger 
diquark-diquark coupling g leads to a stiffer EoS. A diquark-diquark coupling 
g induces enough repulsion in order to prevent the collapse of the system to 
infinite diquark BEC condensate and makes the system stable. The inserted 
plot in Fig. 6 shows more clearly the appearance of several branches of P{pb) 
which is a manifestation of first-order phase transition. At a fixed baryonic 
chemical potential only the highest positive pressure corresponds to a sta- 
ble phase. In Fig. 6 (left pannel), we show the baryonic density pb/po as a 
function of the baryonic chemical potential pi, for various parameter sets. The 
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Table 4 

The parameters Ajy, and 5 in the deconfined phase of the model when the critical 
deconfinement chemical potential = 923.25 MeV is assumed. Parameter set B2 
is already given in table 3. All parameter sets lead to baryonic density = pQ at 
the critical deconfinement chemical potential p^ = 923.25 MeV and reproduce the 
empirical saturation point {Eh/ A = —15.75 MeV, po = 0.15 fm^). 
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Fig. 7. The binding energy per baryon Eb/A as a function of density Pb/ Pq for the 
nuclear matter with parameter set NM (table 1) and for the diquark matter with 
different diquark ultraviolet cutoff A/) given in table 4. 

deconfinement phase transition leads to a baryonic density jump at the critical 
chemical potential ^\ = 923.25 MeV similarly to the nuclear matter liquid- 
gas phase transition. In this case, the deconfinement transition can only be 
first-order since in hadronic side, the pressure has multi-branches. 

In order to understand the implication of different choices for the diquark 
ultraviolet cutoff A^, we also calculated the diquark EoS for three different 
values of the cutoff A^i given in table 4, assuming that the deconfinement 
takes place at /i* = 923.25 MeV. All parameter sets given in table 4 repro- 
duced the empirical saturation point {Eb/A = —15.75 MeV, po = 0.15 fm^). 
It is seen from table 4 that by increasing the cutoff A^), a weaker repulsive 
diquark-diquark coupling g is needed. In Fig. 7, we show the binding en- 
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ergy per baryon for the hadronic matter with parameter set NM (in table 1) 
and for diquark matter with various parameter sets given in table 4. For all 
diquark parameter sets in table 4, at the critical deconfinement chemical po- 
tential fj^l the baryonic density is equal to the normal nuclear matter density 
Pb — Po — 0.15 fm~^. It is obvious that at high density diquark matter is 
energetically favored compared with nuclear matter while at low density nu- 
clear matter is the stable one. Note the EoS of the hadron or the diquark EoS 
alone has a minimum at the saturation density since we have = Pq = 
at nl = 923.25 MeV, reflecting the thermodynamical consistency equation 
(37). The produced nuclear-diquark EoS is softer than the pure nuclear mat- 
ter EoS^. By lowering the diquark ultraviolet cutoff A^, the diquark EoS 
becomes stiff er, see Fig. 7. 

It is plausible that despite the fact that we have turned off the vector-meson 
field at the saturation point, the diquark interaction induces enough repulsion 
to stabilize matter. In other words, in our model diquarks play the role of the 
vector-meson in the deconfined phase and the rearrangement of quarks from 
packages of three quarks to packages of two quarks is essential. In a quark 
exchange model of NN-interaction, the s-channel picture of quark exchange 
looks like two remaining diquark cores and a diquark in between. Therefore, 
one could say that the diquark picture is dual to a meson exchange picture. 
The bosonic approach reproduces the correct saturation point and a soft EoS. 
In fact, repulsive bosons are easier to handle than fermions where the corre- 
lated wave function cannot be adequately described by the mean-field theory. 
Moreover, the diquark BEC goes over directly into BCS-condensate. 

In Fig. 8, we show the energy per baryon eh,q/pb as a function of inverse 
baryonic density {pb/Po)~^ ior two parameter sets given in tables 1 and 3. 
They exhibit a typical Van der Waals-like isotherm diagram where hquid and 
gaseous phases coexist. At low density in the confined region, nuclear matter 
is mechanically unstable since the pressure is negative. But the pressure in 
vacuum and at saturation density po is equal to zero. Therefore, the nuclear 
matter mechanical instability can be prevented by fragmenting nuclear matter 
into droplets of nuclei, where each droplet has pressure zero and density equal 
to Po — 0.15 fm~^. Diquark BEC matter also breaks up into stable droplets 
in which the pressure is zero and density is pf. The diquark BEC droplets are 
therefore surrounded by non-perturbative vacuum and droplets of nucleons. 
As we already discussed at the chemical potential pb = 923.25 MeV chiral 
symmetry is restored. Therefore, diquark droplets exhibit QCD-matter in a 
chirally restored form, whereas the chiral symmetry within nucleonic droplets 



^ A computation of the compression moduhis of the joined miclear-diquark EoS 
(solid and dashed lines in Fig. 7) exactly at the saturation point is not possible, 
since we have assumed a sharp boundary between the baryonic and the diquark 
phases. 
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Fig. 8. Energy per baryon number as functions of inverse baryonic density {pb/po)~^ 
in the unit of the nuclear matter density pQ. Various parameter sets given in tables 
1 and 3 are used. The dotted line represents the Maxwell-construction which deter- 
mines the region of coexistence: between the nucleons gas, hadronic droplets with 
7^ 0, the droplets of diquark Bose-Einstein condensate with a = surrounded by 
a non-trivial vacuum with a = 93 MeV. 

is only partially restored. The nucleon mass at the saturation density is about 
763 MeV for the parameter set NM, see Fig. 4. Therefore, in our scenario, 
multifragmentation has more facets than usually. It also signals a first order 
phase transition of chiral restoration and deconfinement. Ordinary nuclear 
matter can then be identified as being in the mixed phase of the first order 
phase transition: consisting of droplets of diquarks and droplets of nucleons 
surrounded by empty space. 



Within this model we cannot directly identify the diquark BEG droplets with 
abnormal nuclei, since color should be gauged if the model is to yield color- 
singlet droplets. However, one may tempt to do so, since we know that within 
nucleons the chiral symmetry is already restored, and the quark density is 
not zero. In this case, one arrives at a similar picture to that in the MIT bag 
model. A similar picture of a nucleon as a bag within which chiral symmetry 
is restored has also been proposed by Buballa [30], Alford, Rajagopal and 
Wilczek [11] and Berges and Rajagopal [12], (see also Ref. [31]), but within 
quite different models and mechanisms. Note that here we can easily choose 
parameter sets where baryon density in droplets of diquarks is different from 
the nucleonic droplet ^ po- Identifying a diquark droplet with abnormal 
nuclei is possible for p{,i > po? this may favor a soft EoS, see Fig. 6. In this 
picture, it is also possible to have combined droplets of nuclei and diquarks. A 
fraction of baryonic charge can then be delocalized in form of a diquark BEG. 
In this case, diquark concentration within nuclei can be smaller than nuclear 
matter density po. A possible chemistry of such a mixed matter is beyond the 
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Fig. 9. Pressure as functions of Pb/ Po for various parameter sets given in tables 1 
and 3. Parameter sets Al and A2 (table 2) are given with the assumption that the 
deconfinement takes place at p*" = 1214.5 MeV, corresponding to baryonic density 
Phi Po = 3. Parameter set A2 leads to a coexistence region between nuclear and 
chirally restored diquark BEC matter within 3 < Pb/ Po ^ 4. Parameter sets Bl 
and B2 (table 3) are given with the assumption that the deconfinement takes place 
right at the nuclear matter saturation point p^ = 923.25 MeV. The gray closed line 
corresponds to the Danielewicz et al. constraint [29]. 

scope of this paper. 

There are some experimental restrictions coming from the flow analysis in 
heavy-ion collisions which indicates the preferable theoretical values of pres- 
sure in a portion of baryonic density 2 < Pb/po < 4.6 at zero temperature 
[29]. This region is shown in Fig. 9 with a gray closed line. It is important 
to note that no single EoS can simultaneously describe all phenomenological 
data [29]. The area of allowed values does not represent experimental values 
itself but results from transport calculations for the motions of nucleons in a 
collision for a given EoS. To compare with the flow constraint region, we sepa- 
rate regions of SIS [pb < 3po) and AGS {pb > 3po) energies as weak and strong 
flow constraints. In this figure, we show the EoS we obtained for two cases of 
early deconfinement (sets Bl and B2 in table 3) and late deconfinement (sets 
Al and A2). In both cases, the produced equations of state pass through the 
acceptable phenomenological region. For the region pb < 3po the EoS obtained 
by assuming rather early deconfinement (sets Bl and B2 in table 2) passes 
through the constraint region. These equations of state stay outside the ac- 
ceptable phenomenological region at higher density. However, because of the 
diquark ultraviolet cutoff, at higher density of about pb > 2.6po, the model 
is not trustworthy. The current available phenomenological data from heavy 
ion collisions cannot limit the acceptable boundaries on pressure as a function 
of density below two times the normal nuclear matter density [29]. Note that 
the models (for sets Bl and B2) assuming the deconfinement phase transi- 
tion at relatively low baryonic density cannot be precluded by the present 
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flow constraint. This is in contrast with conclusion made in Ref. [29]. For 
the region pf, > 3po the EoS obtained by assuming late deconfinement passes 
through the allowed region. For parameter sets Al and A2, we assumed that 
the deconfinement takes place at baryonic density pi, = 3po- For parameter 
set A2, the deconfinement phase transition is first-order and it leads to a co- 
existence region between nuclear and chirally restored diquark BEC matter 
within 3 < pb/po < 4. For parameter set Al, the deconfinement phase tran- 
sition is second-order and there is no coexistence region. Note that since we 
used rather stiff EoS for the hadronic side (set NM) , in the confined region the 
nuclear matter EoS stays outside allowed region. This is a general problem of 
these simple mean-field models and may be improved by adding further meson 
interactions or potential terms to the Lagrangian [23]. 



6 Summary and conclusions 

In this paper, we investigated the nuclear matter and its possible phase tran- 
sition to diquark BEC matter at zero temperature. For the hadronic phase 
we used a linear a — uj model incorporating explicitly the effect of Dirac-sea. 
For the deconfined phase, we used a complex scalar field theory in terms of 
diquarks which couple to mesonic fields and repel each other. Nucleon, me- 
son and diquark masses are dynamically generated via spontaneous symmetry 
breaking. 

Our simplified theory does not yet give uniquely the order of a possible de- 
confinement transition, nor its location. However, as pointed out, there are 
indications which suggest that it might be first order and takes places above 
nuclear matter saturation density. In this paper, we studied the implication of 
different choices for the critical deconfinement chemical potential. Our knowl- 
edge about the form of matter above deconfinement is still very rudimentary. 
In this paper, we studied a possible formation of a diquark BEC phase with 
deconfinement and investigated its implication for the baryonic matter. Since 
we do not know the strength of the diquark-diquark interaction coupling g, 
we investigated possible values for this coupling and possible physical conse- 
quences. 

We found that it is possible that nuclear matter undergoes a diquark BEC 
phase transition. Our study generally indicates the importance of the diquark 
BEC phenomenon at a rather low chemical potential. It has been shown within 
various models that phase transition from hadronic matter to normal quark 
matter cannot occur [13,33]. Based on our studies, deconfinement to diquark 
BEC matter seems to be a viable option. 

We showed that chiral phase transition is of first order and coincides with 
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deconfincmcnt transition. The compressibility of the EoS depends on the 
diquark-diquark couphng and the critical deconfinement chemical potential. 
We found that an early deconfinement and a weaker diquark-diquark inter- 
action (or a stronger diquark BEC) soften the EoS. Assuming that decon- 
finement already takes place at saturation, we proposed a scenario in which 
the satTiration properties can be described by the diquark BEC phenomenon. 
In contrast to the standard picture, the vector-meson is not needed and the 
diquark-diquark interactions provide enough repulsion in order to prevent the 
collapse of the system. The EoS obtained in this way are rather soft. 

In particular, we studied the various coexistence patterns which may emerge if 
deconfinement takes place at different critical chemical potential. We showed if 
deconfinement takes place at higher baryonic chemical potential, due to a first- 
order phase transition there exists a coexistence region between the hadronic 
and the diquark matter. The extension of this region in density depends on the 
diquark interaction, a weaker diquark-diquark interaction widens this region. 
For model parameters with heavy diquark in vacuum (albeit within baryons), 
phase transition to diquark matter becomes of second-order. We also showed 
that if deconfinement takes place at about the nuclear matter saturation chem- 
ical potential, the ordinary nuclear matter can then be conceived as being in 
the mixed phase of the first order phase transition, consisting of droplets of 
nucleons and droplets of chirally restored BEC diquarks surrounded by a non- 
perturbative vacuum. 

Diquark Bose-Einstein condensation guarantees a substantial presence of di- 
quarks in the plasma. A heavy quark-antiquark pair (QQ) may be converted 
into a heavy baryon (Q-diquark) and exotic unstable pentaquark (Q-diquark- 
diquark) in the presence of the diquark BEC. Also the ratio of Ac/S^. can be 
altered in presence of the diquark BEC, since Ac contains a scalar diquark 
while the production of Ec remains unaffected by the presence of the scalar 
diquarks [34]. It will be challenging for incoming facility at GSI/FAIR to probe 
a possible signature for the diquark BEC matter. Moreover, the diquark BEC 
may alter the compact star cooling behaviour. In principle, two quarks can 
be converted into a diquark and a pair of neutrinos via the following weak 
interaction process u + d ^ [^^o^] diquark + i^e + ^e- These extra neutrinos will 
be strongly enhanced in the presence of the diquark BEC. This mechanism 
accelerates the cooling of compact star. On the other hand, it has been shown 
that the cooling process will be suppressed in the presence of BCS-pairing 
[35]. Therefore, the cooling behaviour of compact star may provide a possible 
signature of existence of the diquark BEC matter and BEC-BCS crossover. 

It is highly desirable to improve our knowledge about the nature of diquark 
interaction in medium and its possible connections to QCD parameters or some 
phenomenological inputs. It is also of interest to investigate the implication of 
the nuclear-diquark model presented here under neutron stars conditions. 
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